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Motivated by examples of erosive incision of channels in sand, we investigate the motion 
of individual grains in a granular bed driven by a laminar fluid to give us new insights 
into the relationship between hydrodynamic stress and surface granular flow. A closed 
cell of rectangular cross-section is partially filled with glass beads and a constant fluid 
flux Q flows through the cell. The refractive indices of the fluid and the glass beads are 
matched and the cell is illuminated with a laser sheet, allowing us to image individual 
beads. The bed erodes to a rest height hr which depends on Q. The Shields threshold 
criterion assumes that the non-dimensional ratio 9 of the viscous stress on the bed to 
the hydrostatic pressure difference across a grain is sufficient to predict the granular 
flux. Furthermore, the Shields criterion states that the granular flux is non-zero only for 
9 > 9c- We flnd that the Shields criterion describes the observed relationship hr cx Q^^^ 
when the bed height is offset by approximately half a grain diameter. Introducing this 
offset in the estimation of 9 yields a collapse of the measured Einstein number q* to a 
power-law function of 6* — 6'c with exponent 1.75 ± 0.25. The dynamics of the bed height 
relaxation are well described by the power law relationship between the granular flux 
and the bed stress. 



1. Introduction 

The response of a granular bed to forcing by a fluid which flows through and over 
the be d has been t he su b ject o f continuous inquiry for over a century (see, for ex- 
ample Graf (1971); YalinI ( 1977? )). This phenomenon is at the centre of a wide range 
of practical and fundamental problems. Predicting the granular flux for a know n fluid 
flow is important in u nderstanding h ow beaches ( Bailard 'l98l'; 'Komar' 'l99^, rivers 
(jMurrav fc Paolalll994l) and delt as (Ke nvon fc Turcot tc 1985; Parker et a?. 19 98,) evolve, 
mountains erode ( Burbank et al\\l99(h and landscapes form (Howard! il99l . Sedimen- 
tary records cannot be decipher ed without a worki ng understanding of the combined 
fluid/granular (two-phase) flow ( Blum fc Torngvis^ ,2000). The microscopic details of 
the bed's response to forcing by a flui d raise impo r tant fundamen tal questions about 
the nature of fluid fl ow near a rough ('Grassl ll97ll: [j imenea 120041 ) and/or permeable 
( Brinkmanl 194 9: Be avers fc Jos eph 1967) wall, the motion of a grain on a rough sur- 
face (ISamson aLlll999HQuartier et aLll2000l) . and the dynamics of granular avalanches 
(jPouadv et a/.l 120021 : iBonamv et al\\20o£ . 

Given some measure of the fluid flow intensity, one would like to predict the granular 
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flux from the properties of the granular material such as size, shape, friction cocfflcient, 
bed packing, etc. An important aspect of the problem to consider first is the onset of 
granular flow. J ust as the surface o f a granular pile driven by gravity alone relaxes to the 
angle of repose ( Jaeger et al. 1989I) . granular beds driven by fluid flow are thought to be 



static below a certain threshold fluid flux. Several empirical curves relating some measure 
of the fluid forcing at the onset of granular fl ow to the grain properties have been pro- 
posed ( Shields! 19361 : Vanoni|[l964l : lYalin 1977). Numerous experimental studies aim ed at 
compu ting these curves (^summarised in lMiller et al. ( 1977 ) or feuffing t on fc Montgomery! 
(|l997l ). for example) are diflicult to interpret since measured critical properties (such as 
the Shields parameter) vary by as much as a factor of three from study to study. The dis- 
parate and subjective definitions of the onset of g ranular fiow used by the resea, r chers are 
frequ ently used to explain the observed scatter (jShvidchenko fc Pendei![2OO0l : IPaphitid 
20011) . 

Even the existence of the threshold fiuid forcing b elow which the gran ular flux is iden - 



tically zero is a subject of a lively ongoing debate (jGraf fc Fa zis 197 71: Ivan Riin 



When the flow is turbulent, the local bed forcing is stochastic ()Papanicolaou et al. 



198S) 



2001) 



and a strong case, supported by data, can be made for the prese nce of a granula r flux 
(albeit vanishingly small) for any mean value of the fluid forcing ( van Riin 1989! ). The 
omnipresen ce of a granular flux poses a two-fold problem. First, an arbitrary t hreshold 
gram flux (iNeill fc Yalin!ll969h. d islodg ing rate (jShvidchenko fc Pendei! lioooh . or dis- 
lodging probabilitv (jPancev et al. 2002!) m ust be introduced to characterise the onset. 
Whereas the study of Graf fc Pazid ( 1977 ) suggests that the granular flux is exponen- 
tially small below some threshold driving, no systematic attempt, to our knowledge, has 
been made to quantify the transition in the functional dependence of the granular flux 
on driving. Second, the ever-present granular flux results in the evolution of the sur- 
face packing. This "ageing" or "armouring" of the granular bed leads to a decrease in 
the granular flux and an increase in the perceived thresho ld for the onset of persistent 
granular flow (jCharru et a/.!l2004!: IPaphitis fc CollindbOOSi ). 

Even when the flow is laminar, the rough surface of the granular bed yields a fluctuating 
local fluid stress, albeit deterministically related to the time dependent realization of the 
surface packing. Statistical methods are therefore necessary to completely characterise 
the onset of granular flow. Because we expect the surface packing to have a transient 
component, so will the granular flux. However, we expect there to be a true threshold 
forcing below which there exist surface packings such that the associated surface stress 
is insufficient to dislodge any grains. Under steady pre-threshold forcing conditions, one 
such surface packing will eventually be realized and the granular flow will cease. We are 
interested in predicting the value of the threshold fluid forcing, the transient granular 
flux in pre-threshold conditions, and the steady state granular flux above the threshold. 
We would also like to establish whether a single characteristic of the fluid forcing (such 
as the surface stress) is sufficient to predict the steady state granular fiux. 

Here we report our initial investigation of these questions in a cell partially filled with 
transparent spherical glass beads driven by an index-matched fluid. The small size of the 
cell and the grains and the high viscosity of the fluid ensures that the particle Reynolds 
number is always less than ten. Adding small fluorescing tracers and illuminating the 
cell with a laser sheet allows us to directly image all grains in a vertical slice through 
the system. We can therefore simultaneously measure the fiuid and granular fluxes, the 
dimensions of the fluid filled region, and the granular packing. This allows us to test 
theoretical predictions of the threshold condition and the functional dependence of the 
granular flux on fluid driving. Our study may be considered complementary to that of 
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Figure 1. (a) Schematic diagram of the experimental apparatus, (b) Sample image of the refrac- 
tive index matched beads in a plane away from the side walls, (c) Side view of the experimental 
apparatus. Solid line (1) is the initial profile of the granular bed after the preparation steps 
(i)-(iii) are carried out. Dashed line (2) is a sketch of typical profile of the granular bed after 
the cessation of erosion, the profile is perfectly fiat except near the ends. The image region is 
located about 90 mm (130d) from the upstream end. 



Goharzadeh et al\ (|2005r ). who examined the fluid flow near a granular bed using a similar 
setup. 



2. Experimental technique 

As shown in figure [U an enclosed cell with a rectangular cross-section is partially filled 
with glass beads with diameter d = 0.70 ± 0.05 mm, density pg — 2.5 g/cm'^ and index 
of refraction of n = 1.54. A peristaltic pump is used to generate a constant fiuid fiux Q 
through the cell. The fluid is manufactured to have the same refractive index as the glass 
beads and a dye is added which fluoresces with a peak intens ity at a wavelength of 565 



nm when illuminated by light with a wavelength of 532 nm ( Siavoshi et al. 2006[ ). The 
hydrocarbon based fluid, manufactured by Cargille, Inc, has density Pf = 1.026 g/cm^ 
and kinematic viscosity — 0.243 cm^/sec. 

The granular bed is prepared using the following three steps, (i) we tilt the cell ver- 
tically with the fluid input end at the bottom, allow the grains to come to rest and (ii) 
slowly set it down horizontally. This procedure fills the cell top leaving an empty buffer 
space near the exit mesh (see figure [TJ;). Step (iii) of bed preparation involves applying a 
small fluid flux to create a small (3-4 grain diameters) uniform gap. The three steps above 
are repeated to create the initial condition for all runs. We scanned the laser sheet across 
the initial bed surface and verified that it was fiat in the observation window of width to 
within half a bead diameter. After the bed is prepared, the fluid flux is instantaneously 
raised to some value Q and kept constant for the duration of each run. After traversing 
the pile's surface, eroded grains fall out of the way into the empty buffer space between 
the pile and the downstream mesh. Eventually, all granular motion ceases. In such a rest 
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state the profile is perfectly flat except in regions of order 2 cm near the inlet and outlet 
meshes. We therefore suppose, and argue further below, that in the flat region, which 
comprises at least a 20 cm portion of the cell, the fluid flow proflle is steady, laminar and 
independent of the downstream distance. We further checked that systematically varying 
the initialization fluid flux to create different initial bed heights did not change the final 
bed height. We observed that the final height reached was the same to within a grain 
diameter, provided the initial bed height was greater. 

The range of fiuid fluxes Q for which the rest state is reached is a function of the cell 
filling fraction. If the cell filling fraction is high, the range of Q is rather small. However, 
if the filling fraction is small, the quiescent bed does not have a flat region when the rest 
state is reached. We therefore set the uniform bed level at roughly 22 mm. 

A laser and cylindrical lens system, placed above the cell, illuminates a vertical slice 
through the cell away from the sidewalls. Because the beads do not contain the dye, 
they appear dark against a bright background fluid. A high speed digital camera with 
a resolution of 512x480 pixels records a sequence of images at a rate of 30-60 frames 
per second. A typical image of a vertical slice through the bed is shown in figure [Dd. 
The image is truncated at a depth of about 15 mm (measured from the top plate) below 
which the particles never move. The apparent size of a bead depends on the distance of 
its centre from the illuminated plane. We are able to identify all beads whose centres 
are within approximately 0.3d from the illuminated plane. The height of the granular 
bed in every image is measured using an edge detection algorithm. We take the height 
of the granular bed to be the infiection point in a depth profile of the pixel intensity 
averaged over the image width. This technique yields the mean height of the bed in the 
observation window. Furthermore, the mean height is averaged over 5 experimental runs 
to reduce statistical fluctuations, and remove any variability in preparing the initial pile. 

We verified that the bed height decreases monotonically away from the side walls up 
to a distance of about 3.5 mm beyond which it remains constant to within the roughness 
scale of the granular bed. For all subsequent measurements we therefore positioned the 
laser sheet at a distance 7 mm from the side walls. The imaged region is about 16 mm 
wide, and is located 90 mm downstream from the inlet mesh (see figure [Tt). Comparing 
the height data for experimental runs repeated under same preparation conditions, we 
observed the standard deviations to be less that one particle diameter. 

The inlet and outlet ramps ensure a smooth fiow of the fiuid onto the granular bed. 
After a process of trial and error, the height of the ramps (22 mm) was chosen to be 
slightly lower than the height of the granular bed when the grains fill the cell uniformly 
forming a flat surface. The leading edge of the bed tends to align with the top of the 
inlet ramp. If the bed is eroded below the upper edge of the inlet ramp, an undesirable 
turbulent region can form behind it. Neglecting end effects is theoretically justified be- 
cause we find a relationship between the grain fiux and the local viscous shear stress in 
a region away from both ramps. In addition, we find that the final rest height of the bed 
away from the ramps depends only on ffuid flux Q and not on the bed preparation or 
on the geometry of the ramps. Therefore, we maintain that the ramps only affect the 
granular bed in their vicinity and do not influence the bed height further than a certain 
distance ib downstream. To estimate ib, we identify the length scale at which the width 
6{x) of the Blasius boundary layer equals half the fluid depth hf. The characteristic flow 
velocity is J7 ~ Q/Whf, where W is the width of the cell. Therefore 

4^^^^. (2.1) 
For the largest values of the flux Q and the fluid depth hf in our experiment, the length 
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Figure 2. The horizontal granular velocity (vx) (bottom axis, circles), computed in the regime 
when grains are rolling on the surface of a quiescent pile, averaged over the observation window 
and a one second time interval. The granular volume fraction (j) is plotted with diamonds (top 
axis). The vertical axis denotes the depth z measured from the cell's top plate; z — corre- 
sponds to the top plate; the bed's surface is located where (f> vanishes. The total granular flux 
is proportional to the z integral of the product of {vx)4>- 



£b ~ lOOrf. Our use of the laminar flow profile in the window of observation located ISOd 
downstream from the inlet ramp is therefore justified. 

To compute the vertical profile of the particle velocity and packing fraction we measure 
the position of every particle in the image to within O.ld using a centroid technique 



tne position or every particle in tne image to witmn u.ia using a centroia tecnnique 
( Siavoshi et all [20061 Tsai et al. 2003[) . The velocity Vx of individual particles is then 



determined by tracking the particles over two successive images. The image region is 
divided into horizontal bins of height Id. The volume fraction (f> is determined by dividing 
the number of particles in each bin by the volume of the bin (Id x 0.6d x 22d). The 
depth dependence of the bin-averaged and time averaged (over approximately one second) 
velocity (vx) and volume fraction (j) for one particular flow rate is shown in flgure[2l Above 
a certain depth no grains are detected in the bin which yields a null packing fraction. 
Given the depth profiles of the horizontal velocity and volume fraction, the grain flux qg 
(averaged over the observation window) is obtained by computing the integral / {Vx)4> dz. 
Every reported flux value represents the average over a time interval of one second within 
which the variation in the flux values is small. All the experiments were repeated five 
times and the reported values are the averages over these five runs. 

Before we launch into the detailed analysis of the data, let us qualitatively describe the 
phenomenon. Immediately after the fiuid flow is switched on, the granular bed experiences 
a brief period (a few seconds) of rapid shear flow during which several grain layers are 
moving. The fluid gap increases rapidly during the shear flow regime. Subsequently, the 
granular flux is conveyed via "bed-load" , i.e. individual grains rolling on the surface of 
an apparently quiescent bed. After several minutes the granular flux eventually vanishes. 
Because the erosion is gradual, the fluid gap increases gradually. Thus, the bead shear 
stress decreases gradually and the bed (away from the ramps) must approach a trivial 
flat state in which the fluid exerts a threshold shear stress on the bed. 
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Figure 3. (a) The average depth ft/ of the fluid layer above the granular bed. Solid lines are a 
guide to the eye. (b) The rest depth hr of the fluid layer measured after cessation of granular 
flow as a function of fluid flux Q. Dashed line is a flt to equation (|3.2|l . Solid line is the same 
flt with hr offset by roughly 0.4 mm. For reference, the fluid depth above the entrance ramp is 
11mm. 



3. Rest height of the granular bed 

Figure [3^ shows the average depth of the fluid gap as a function of time for a constant 
fluid flux Q over the bed. As the grains are eroded from the bed's surface and deposited 
into the buffer space, the fluid gap grows, and the mean fluid flow speed and therefore 
the shear stress on the bed decreases. The erosion eventually ceases and the fluid gap 
approaches a constant rest depth which decreases with the applied fluid flux Q. 

The graduaUy increasing fluid gap and thus gradually decreasing bed stress leads to 
the cessation of erosion precisely at the threshold driving. The variation of the rest depth 
hr with the appUed fluid flux Q allows us to examine the validity of the Shields threshold 
criterion. We estimate the bed shear stress a assuming two-dimensional laminar flow 
between infinite smooth parallel plates separated by ft./. This approximation is good 
when the cell width « 26 mm is much greater than the fluid gap depth hf. Using 
the solu t ion fo r laminar flow in a pipe of rectangular cross-section (see for example 
Cornish (fl928h ) we can compute the correction to the estimated bed stress. Neglecting 



the sidewalls results in an underestimate of the stress by 8% for the smallest fluid gap and 
by 29% for the largest fluid gap in our experiment. However, neglecting the sidewalls has 
the advantage of a simple analytic expression for the bed stress which greatly simplifies 
the analysis. The laminar flow assumption is good since the cell Reynolds number for 
Q = 1000 cm'^/min is Re = v^es-nhf/v = Q/Wv « 0.3, where Wmoan is the horizontal 
fluid velocity averaged over the fluid gap. The Shields parameter 9, defined as the stress 
a scaled by the hydrostatic pressure difference [pg — Pf)gd across the grain, is 

(3.1) 



W-/gdhj '' 



where 7 = Pg/pf — 1 « 1.437 is the density contrast. The Shields parameter measures the 
relative importance of the destabilising hydrodynamic forces and the stabilising gravity. 
Assuming that the rest depth hr of the ffuid gap corresponds to the Shields criterion 
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Figure 4. Einstein number q* vs. the excess 9 — 6c of the dimensionless shear stress. Solid line 
is a power-law fit to the region 9 — 9c < 0.3. Dashed line is a power law fit with A = 1.6 fixed. 



9 = 9c, we obtain 

hr = d\ — ] , with Qr ^ 9c — . (3.2) 

\Qr J &V 

The fit to the measured rest depth (shown as a dashed line in figure ^p) is signifi- 
cantly improved if a positive constant of approximately half a diameter d is introduced 
to the right-hand side of l|3.2p (an even better fit is obtained if K ~ Q°-^^). The fit 
yields Qr — 15.4 ± 0.3 cm^/min which translates to the critical Shields parameter of 
9c = 0-30 lb O-Ol- Thi s va lue is consistent with previ o usly r eported values, for example, 
m iMiller et al\ (|l977l ) or iBuffington fc Montgomery! (|l997( ). For comparison, the Yalin 
parameter, S = (^gd^Y/"^ /v, a less widely used dimensionless group which measures 
the relative importance of viscous and gravitational forces without reference to the flow 
intensity, is S « 0.03 in our experiment. 

The physical origin of the offset is not entirely clear to us at this time. The correction 
to the bed stress estimate due to the sidewalls cannot be the source of the offset since 
in the limit of small gap /i/ the sidewall effect vanishes. A possible source of the offset 
is the fluid boundary condition at a permeable wall which relates the derivative of the 
fluid velocity to its magnitude at some penetration length scale t. Yet another possibility 
is that the offset could be due to the way in which /i/ is measured. We use the mean 
surface height to represent the mean boundary stress. In reality, however, fluid flow past 
a rough permeable medium creates a distribution of bed stresses and the estimate of its 
mean using a non-slip boundary condition at a flat wall is a gross approximation. The 
fact that the error in the estimate of the mean stress translates to a height offset that is 
only half a bead diameter is a pleasant surprise. 



4. Flow rule 

Tracking all grains in a vertical slice as mentioned in section [2] allows us to directly 
measure the grain flux qg and to attempt to relate it to the estimated mean bed stress. 
In general, other factors besides bed stress, such as the flow history, could influence the 
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grain flux. Although an unambiguous relationship between grain flux and bed stress is 
usually tacitly assumed, here we test it directly. 

Following convention ([Einstein Il950l) we use the Einstein number q* — qg/{jgd'^Y''^, 



i.e. the non-dimensionalised grain flux. Figure [3] is a plot of q* vs. the estimated excess 
bed stress — Oc computed using equation (|3.ip with the value of Oc = 0.3 measured 
in Section [31 Different symbols correspond to different fluid fluxes Q. Because the same 
bed stress occurs at different fluid gap depths for different Q and therefore entirely 
different times in the history of the erosion process, data collapse suggests that only 
the instantaneous bed stress (and not its history, for example) is required to predict the 
instantaneous grain flux. Therefore, upon changes (in time or space) of the bed stress the 
grain flux quickly relaxes to its corresponding value. We hypothesise from dimensional 
considerations that the relaxation time scale is d^/i/ « 0.02 sec, which is small indeed. 

The data in Figure [4| do not collapse for high stresses, i.e. early times, because small 
fluid gaps are associated with changes in the mechanisms of granular transport. At early 
times the granular flux occurs primarily as bulk shear of the granular pile, whereas later 
in the history of bed evolution the grain flux is conveyed through isolated grains rolling 
over the otherwise immobile granular bed. The latter form of transport is known as the 
"bed-load" regime. A power law flt to the data collapse yields 

q* K. A{e -Oc)^, with A = 1.89 ±0.25, ^ = 0.11 ±0.03. (4.1) 



Our m easured exponent A is slightly higher than t he widely accepted valu e of 1.5 ([Mever-Peter fc Miiller 



19481) . and also above the value of 1.6 found by [Wong fc PaAe3 (|2006l ) in their reanal- 



ysis of Mever-Peter fc Miilleil ( 19481 ). The p re- fact or A, however, is more than an order 



of magnitude smaller than that quoted bv [Wong fc Parker ( 2006[ ) for turbulent flow. 



The large discrepancy suggests that although the granular flux in the laminar fluid flow 
regime is described by the same functional dependence on the Shields parameter as in 
the turbulent regime, the pre-factor is not. The turbulent fluid appears to be signiflcantly 
more efficient at entrainment and transport of surface grains. This is perhaps not sur- 
prising since for the same mean bed stress, the turbulent bursts provide a more effective 
dislodging process than the smoothly varying bed stress in the laminar flow regime. 



5. Approach to rest height 

The experimentally verifled relationship between the bed stress and the granular flux 
allows us to compute the entire history of the fluid gap depth hf(x, t) (x is the coordinate 
in the flow direction) as it approaches the rest depth hr- The mean thickness of the flowing 
grain layer is negligible compared to hj in the "bed- load" regime. Therefore, we write 
the conservation of sediment as 

dhf _dqg 

'df-^- ^^-^^ 
To study the approach of the proflle to the rest height hr, we introduce ^ = hf/hr, scale 
all lengths by d and times by yd/Tg and use the flow rule (|4.1|) which can be expressed 
as q* — A0^{1/^'^ ^ l)'^- Eq. ()5.1|1 can therefore be written as 

§+c(e)§ = 0, ci0^2XAe^r'-'\l~e)'-'^^- (5.2) 

Given an initial condition ^(x,to) — ^0(2;), the solution to equation (|5.2p (obtained by 
the method of characteristics) determines ^{x,t) implicitly through ^{x,t) — ^o{x — 

ciax,t))t). 

There are two issues which make the full solution of (|5.2p difficult to use. First, a 
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possibility exists that the granular flux includes a term proportional to the gradient of 
the bed height resulting in a diffusive term in ()5.2|) . Our measurements of the bed height 
are not sufficiently resolved to verify or reject such a diffusive term. Second, the domain 
of applicability of equation (|5.2p is limited in time and space. In the beginning of the 
experiment, granular transport is mainly conveyed via bulk shear and thus does not 
obey (|4.ip . Also, near the inlet ramp the bed height may change rapidly and thus the 
approximate expression for the shear stress (j3.ip is invalid there. Fixing the initial and 
boundary conditions for equation (|5.2p is therefore not feasible. Thus the full solution is 
not useful. 

However, the long time asymptotic solution to (j5.2p with a small diffusive term may be 
insensitive to the initial and boundary conditions. Therefore, we seek a separable solution 
to (|5.2p in the limit ^ ^ 1 (i.e. nearly flat bed) 

X — Xq ^ 



where xq and to are integration constants and B = \A{2dc)^d/hr. 

Because the height of the bed changes by less than a diameter along the cell, an 
experimental test of the spatial variation in Eq. (|5.3p is impractical in our setup. We 
therefore focus on the approach of the fluid depth averaged over the observation window 
to the rest depth hr (cf. figure [3^). Equation (|5.3p implies that once the fluid gap depth 
hf is scaled by rest depth hr and time is scaled by hr/^/^/gd, the data should collapse 
onto a master curve which approaches unity as the power law (t — tQY^^^^^\ Figure [5] 
shows the data collapse. The power law fit to the resulting data cloud yields the offset 
time to w —2.2 hr/^/jgd and the exponent A = 1.6 ± 0.1. The predicted data collapse 
suggests that the assumed functional relationship between the instantaneous bed stress 
and granular flux is indeed correct. It is difficult to interpret the value of to, since the 
asymptotic solution (|5.3p only applies at long times. The fitted value for A falls just 
outside t he la interval o f the value quoted in (|4.ip . but is consistent with the value 
found bv lWong fc Parked (|2006h . 



6. Discussion 

In summary, we have revisited the well-studied problem of a sand bed driven by an 
overlying viscous fluid. We focused on precise measurements in the laminar flow regime 
where little data is available. Our flow geometry, in which the depth of the fluid gap 
gradually increases as the grains are eroded, results in the cessation of granular flow 
and thus allows for an unambiguous definition of the threshold condition. Therefore, 
by varying the fluid flux Q through the cell, we have verified the hypothesis that the 
boundary viscous shear stress is sufficient to predict cessation of granular fiow over the 
surface of a horizontal granular surface. The directly measured value of the critical Shields 
parameter 9c — 0.30 ± 0.01 and the Yalin parameter S = 0.03 is consistent with the 
extrapolation of the previously reported values of 9c to this value of S. The advantage of 
approaching the threshold from above is that there is no ambiguity or subjectivity in its 
definition. On the other hand, a relatively rapid approach to the quiescent state in our 
setup precludes a thorough study of the "armouring" phenomenon. Finally, we remark 
that the threshold for the cessation of granular flow may differ from that for the onset 
of granular flow. This effect is difficult to quantify in general since the onset of granular 
flow is marked by transient flow and steady driving is impossible in our cell since any 
erosion leads to a decrease in driving. 
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= 500 cm'^/min 
Q = 750 cm'^/min 
Q = 1000 cm^/min 
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Q = 1500 cm'^/min 
= 1750 cm'^/min 
Q = 2000 cmVmin 
■Fit with A = 1.6 
Fit with A = 1.9 



Figure 5. Fluid gap depth scaled by the rest depth hr vs. time scaled by hr/^/'jgd. Solid line is 
a fit to a function that approaches unity as {t — to)^^'^~^' with A = 1.6 and to = —2.2 hr/^^gd. 
For comparison, we include a fit to a power law with A = 1.9 (dashed line). 



We directly verified that the boundary shear stress determines not only the onset of 
granular flow, but also the granular flux over a roughly horizontal bed in non-steady 
driving conditions. This is an important result since it implies that the granular flux 
quickly adjusts to changes in the driving. We have established the power law relationship 
between the grain flux and the excess boundary shear stress via two independent mea- 
surements. We measure the granular flux directly by particle tracking and, alternatively, 
we analyse the approach of the bed height to the quiescent state. The average of the 
two independently measured exponents is A = 1.75 ± 0.25. The pre-factor in the power 
law relationship between the Einstein number and the excess Shields stress is more than 
an o rder of magnitude sm aller than that in the widely used "bed-load" granular flux 
rule ( Wong fc Parker 20061 ). indicating perhaps a different mechanism for the granular 



transport in the laminar fluid flow regime. 

We have measured the parameters in the granular flux rule (j4.1l) not only for a par- 
ticular range of Reynolds numbers, but also for a particular type of grains — smooth 
spheroids. How the parameters in (14. ip depend on the grain properties such as shape, 
friction constant and roughness is an interesting and open question. 

Lastly, we note that since the threshold condition in our experiment is reached gradu- 
ally, our problem is similar to the relaxation of a dry granular pile expressed in terms of 
the angle of repose. Moreover, the concept of a Shields threshold, typically defined via 
the onset of granular flow due to gradually increased fluid driving, may be analogous to 
the notion that dry granular piles can be characterised by a maximum angle of stability. 

The work was funded by the Department of Energy grants DE-FG0202ER15367 (Clark), 
DE-FG0299ER15004 (MIT), and the National Science Foundation grant number CTS- 
0334587 (Clark). 
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